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Section  I.  INTRODUCTION 


The  cone  is  a  very  basic  aerodynamic  shape,  yet  no  single  theory 
exists  which  predicts  accurately  the  pressure  distribution  on  the  sur¬ 
face  of  a  cone  for  all  Mach  numbers.  Subsonic  compressible  flow 
theories  [1,2,3]  exist  which  predict  accurately  the  pressure  distribu¬ 
tion  on  a  body  of  revolution  for  all  subsonic  Mach  numbers  not  too 
close  to  one.  Supersonic  flow  theories  [3,4]  exist  which  predict 
accurately  the  pressure  distributions  on  a  body  of  revolution  for  a 
wide  range  of  supersonic  Mach  numbers  above  the  shock  attachment  Mach 
number.  Several  approximate  transonic  flow  theories  exist  (such  as 
Yoshihara  [5],  Leiter  and  Oswatitsch  [6],  Spreiter  and  Alksne  [7],  and 
Wu  and  Aoyama  [8]),  which  predict  the  pressure  distributions  on  a  body 
of  revolution  for  sonic  frecstream  flow,  except  Wu  and  Aoyama' s  method 
was  intended  for  various  ranges  of  transonic  Mach  numbers.  A  comparison 
'  of  various  transonic  theories  was  given  by  Wu  and  Aoyama  [9]. 

It  was  observed  that  calculations  of  the  pressure  distribution  on 
a  cone  using  the  previously  developed  method  of  Wu  and  Aoyama  [8]  were  in 
very  good  agreement  with  experimental  data  except  near  the  shoulder  of 
the  cone.  Thus  the  method  is  very  promising  to  extended  to  low  Mach 
number  flow.  Further  investigations  of  that  method  showed  that  at  a 
particular  fixed  location  on  the  cone  surface,  the  variation  of  the 
pressure  coefficient  with  Mach  number  showed  a  correct  trend  for  all 
Mach  numbers  up  to  the  shock  attachment  Mach  number.  The  only  defi- 
ciences  seemed  to  be  at  low  subsonic  Mach  numbers  where  the  experi¬ 
mentally  determined  pressure  coefficients  were  constantly  lower  than 
the  Wu-Aoyama  calculated  values,  and  near  the  cone  shoulder,  where  the 
experimental  values  were  again  lower  than  the  theoretical  values.  Both 
of  these  deficiencies,  however,  were  attributed  to  the  fact  that  Wu 
and  Aoyama  had  neglected  one  term  in  their  solution  because  of  problems 
encountered  in  their  numerical  iteration  scheme. 

The  subsonic  theory,^  on  the  other  hand,  showed  very  good  agreement 
with  data  for  low  Mach  numbers,  but  the  theoretical  pressure  coefficient 
at  a  fixed  position  on  the  cone  surface  displayed  an  incorrect  trend 
with  Mach  number  near  a  Mach  number  of  unity.  This  deficiency  was 
attributed  to  the  Prandtl-Glauert  form  of  the  governing  equation  which 
becomes  inadequate  near  a  Mach  number  of  one. 

In  view  of  the  previous  observations,  it  was  decided  that  accounting 
for  the  neglected  term  in  the  previous  Wu  and  Aoyama  calculation  might 
improve  the  solution  near  the  cone  shoulder  and  perhaps  improve  the  sub¬ 
sonic  solution  as  well.  Therefore,  a  modified  approach  was  taken  for 
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this  investigation.  This  improved  and  more  consistent  lnvlscid  flow 
ttv  r ry  has  been  employed  to  predict  the  surface  pressure  distribution 
on  a  conical  forebody  at  zero  and  small  angles  of  yaw  for  subsonic  and 
transonic  Mach  numbers  up  to  the  supersonic  shock  attachment  Mach 
number.  Unlike  the  original  solution  of  Hu  and  Aoyama,  no  terms  in  the 
present  analysis  were  explicitly  neglected.  In  view  of  the  observation 
of  the  fact  that  the  solutions  agree  very  well  with  data  even  to  a  very 
low  subsonic  flow,  the  current  method,  though  it  is  only  a  first-order 
approximation,  provides  a  means  of  solution  from  subsonic  to  transonic 
Mach  numbers.  Thus,  Incorporating  this  supersonic  conical  theory  when 
the  shock  attached  to  the  body,  one  can  solve  the  entire  Mach  number 
range  with  continuity. 


Section  II.  FIRST  ORDER  SOLUTION  FOR  FLOW  OVER  A  CONE 


The  nondimens Iona 1  transonic  small  perturbation  equation  can  be 
written  as  [ 8] : 


where  all  lengths  are  normalized  with  respect  to  the  axial  length  of 
the  cone  (Figure  1)  and  all  perturbation  velocities  are  normalized  with 
respect  to  the  free-stream  velocity,  U .  Thus, 


u  v  e  w 

$x  “  U  *  *r  U  ’  -  ~  U 

00  00  r  00 


In  transonic  flow  over  a  cone,  the  sonic  velocity  is  fixed  very 
close  to  the  shoulder  and  the  flow  on  the  surface  over  the  entire  conical 
forebody  is  subsonic,  provided  the  shock  wave  is  not  yet  attached  to  the 
cone  [1],  Hence,  the  coefficient  of  can  be  expected  to  remain 

greater  than  zero. 

Using  the  same  linearization  technique  as  Wu  and  Aoyama  originally 
proposed  [8],  let 

f2(x)  =  (l  -  M2)  -  (7  +  1)M2  *-  (2) 

and  assume  f(x)  is  a  known  function.  The  value  will  be  determined  later. 
By  this  process,  Equation  (1)  is  linearized  to  the  form 


2  11 
f  $>--  +  +  — -  $  =  0 

xx  rr  -  r  -2  66 


★The  factor,  p  =  71  •  M  ,  can  be  absorbed  in  the  transformation  of 

coordinates  as  In  the  original  report  [1].  However,  it  behaves  as  a 
constant  in  transformation;  thus,  it  will  not  affect  the  final  result. 


where  K(c)  is  a  higher  order  small  function.  Assume  that  the  complete 
solution  may  be  expressed  as 

♦  =  *o  +  *i  +  •  •  •  *  (6) 

where  is  the  first-order  solution  corresponding  to  K(c)  =  0,  and 
is  the  second-order  solution.  and  satisfy  the  following  equations: 


(7a) 


and 


1 

r  dr 


K(€)4»q 

i 


(7b) 


A  method  of  successive  approximation  can  be  set  up  to  obtain  the 
complete  solution.  In  this  study,  only  the  first  approximation  is 
treated  and  yields  a  very  satisfactory  result.*  From  here  on,  the  sub¬ 
script  0  for  the  first-order  solution  is  dropped  for  simplicity. 

The  contribution  due  to  the  cross  flow  is  assumed  to  be  small  if 
the  angle  of  attack  is  small.  Let  be  the  axial  flow  component  and 

$  be  the  cross  flow  contribution;  then  it  is  assumed  that  the  non- 
c 

dimensional  perturbation  velocity  potential,  <t> ,  is  composed  of 

*(6,r,0)  =  *a(6,r)  +  *cU,?,e)  .  (8) 


The  linearized  cross  flow  solution  for  slender  bodies  of  revolution 
according  to  [1]  is: 


<,c(l»r*6)  -  8*n  a  COB  ® 


(9) 


*It  was  observed  that  the  transformation  function,  f,  is  'early  constant, 
which  accounts  for  a  small  value  of  f';  hence,  good  agreement  with  data 
results  even  for  the  first  approximation. 
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where 


R«>.sm 


and  the  general  solution  of  Equation  (7a)  for  the  axial  flow  is  [  1] : 

F* 


mi 


3  /  2  2^2* 

0  /  (F  -  tr  +  pr 


dt 


(10) 


where  |*  is  the  value  attained  by  5  at  the  shoulder  of  the  cone  (at  x 
1),  and  F(t)  is  the  source  distribution  which  can  be  determined  from 
the  boundary  conditions . 

Van  Dyke  [3]  showsthat  Equation  (10)  can  be  expressed  by  the 
following  first-order  approximation: 


*a(§’?)  =  tn  ?  +  G(l) 


(11) 


where 


C(5)  -  ^  l 


i* 


n— _ r  +  £  [  dt  .(12) 

2^1(1*  “  I)  J  i  *  t| 


Hence, 


0  (F.r)  =  in  r  +  in  — ■■  ^ . — 

a  2"  *  2  j  ,<5* .  r> 


+ 


4n 


f  dt 

o 


(13) 


The  boundary  condition  for  inviscid  flow  requires  tangential  flow 
on  the  surface  of  the  cone.  In  terms  of  the  untransfrrmed  variables, 
this  boundary  condition  is  expressed  as: 


d$a(x,r) 


dr 


on  the  body 


=  <dUx) 

dx 


=  tan  5"  B 


(14) 


9 
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for  small  angles  where  6  is  the  half  angle  of  the  cone  measured  In 
radians . 

Because  r  was  affected  little  by  the  transformation,  the  boundary 
condition  in  terms  of  the  transformed  variables  [10]  becomes: 

^a 

— T  <5,0  *  5  .  (15) 

on  the  body 


Now  the  source  distribution  F(g)  can  be  found  by  differentiation  of 
Equation  (11)  with  respect  to  r. 


d*a(g,r)  _  F(6) 
dr  ,  2nr 

hence, 

_ 

F(5)  =  2rtr  — - 
br 


(16) 


(17) 


Now  on  the  body 


and  from  Equation  (14) 


m  m 

r  =  R 


(18) 


(19) 


which  implies 


R(5)  =  5x(5) 


(20) 


By  substituting  Equations  (15),  (18),  and  (20)  into  Equation  (17), 
we  obtain: 


F(5)  =  2«62  x(§) 


(21) 
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Now  by  substituting  Equation  (4)  into  Equation  (21),  we  obtain: 


F(t)  =  f  f(t)dt 


Substitution  of  Equation  (22)  into  Equation  (13)  gives  the  follow¬ 
ing  expression  for  the  axial  flow  velocity  potential  in  terms  of  the 
unknown  function,  f  ( | )  : 


♦.<«.*>  - SiR^ 


5  2  5 

f  f(t)dt  +  ^  f  f(t)dt  fin 

0  0  L 


2/  1(1*  -  0 


+  ^  f  0 

2p  J 


5*  I  t 

r  f  f(t)dt  -  /  f(s)ds 
0  0 


I  -  t 


Combine  Equations  (20),  (4),  and  (9);  the.i  the  cross  flow  velocity 
potential  may  be  expressed  in  terms  of  f(|)  by: 


0  (t,r,6>)  -  sin  a  cos  6  yr  f  f(t)dt 

5  0 


Combination  of  Equations  (8),  (23),  and  (24)  provides  an  expression 
for  the  total  velocity  potential  which  satisfies  Equation  (7a). 


*U,r,e)  = 


(in  r)62 


f  f(t)dt  +  ^-  f  f(t)dt  Tin 

0  0 


2/T(i*  -  i) 


L  f  o 
23  J 


5*  e  5 

r  /  f (t)dt  -  J  f (s)ds 

/  0  0 


i  -  t 


2  ^ 

+  sin  a  cos  6  -y-  f  f(t)dt 

Lo 


•I 


«.  .-A 


•'  —  A 


Now  one  is  in  a  position  to  solve  for  the  unknown  function,  f(£)>  Froir 
Equation  (2), 


*x 
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(26) 


and  from  Equation  (4), 


ii  _ 

dx  '  f<*> 

Thus,  the  chain  rule  gives: 

do  _  _d£ 

dx  ^  dx 


or 


(27) 


(28) 


do 


jX 


6  do 
f(!)  oi 


(29) 


Another  expression  for  o-  can  now  be  obtained  by  differentiation 

of  Equation  (25)  with  respect  to  x.  Then  by  equating  the  two  expressions 
for  0-  an  integral  equation  for  f(£)  is  formed. 

For  the  sake  of  orderliness  when  differentiating.  Equation  (25) 
can  be  rewritten  as  follows: 


*(6»r»0)  =  $a  +  o  +  o  +  o  , 

al  a2  a3  c 


(30) 


and  then  one  can  differentiate  each  term  separately  to  obtain: 


x  f(|) 


&a  &a  ^  1 
al  a2  a3  ^c  I 

“ST  +  ~  +  ~  +  ~h\ 


(31) 
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Term  $  : 

_  3 


l* 


:A  / 


fi  f(t)dt  -  J5  f(a)ds 


1 5  -  t| 


dt 


(38) 


or 


Ui)  it* 

"  20  i0  ’ 


(39) 


where  1^*  is  the  subsonic  slender  body  integral  defined  by  Van  Dyke  [3]: 


if  =  f*  nsLussi  dt  , 
0  I  le-t 


and  the  derivative  of  this  Integral  is 


(40) 


■  !*' 


/ 


dt  +  f(£) 

u-t|  5“° 


4  *0?)  -  r.ca 

+  6*  -  | 


(41) 


Now, 


c*. 


Si  _  L.  jt*' 
"  20  0 


(42) 


From  Equations  (4),  (38),  (39),  and  (40), 


i 


F(0  =  f  f(t)dt 


t* 


F(5*)  =  f  f(t)dn  =  0 
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Consider  now  the  terras  on  the  right-hand  side  of  Equation  (30): 

Term  4)  : 

_ *1 


(In  r)52 

3 


I 

/  f(t)dt 
0 


differentiation  gives: 


81  _  (in  r)52 

t  ~  P 


HO 


Term  4>  : 

a2 


now. 


and 


therfore, 


—  P  1 
2/ 6(8*  -  0 J 


=  -Jl*-  2» 
2 id*  -  o 
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(32) 


(33) 


(34) 


(35) 


(36) 


*2 

Hi)  jin 

S2(i*  - 

2PSCS*  -  {) 


6 

2  J  8(5*  -  0 
i 

j  f(t)dt 
0 


(37) 
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F'(i)  =  f(0 


F (t)  =  J  f(»)ds  , 

0 

F'(t)  =  f(t) 

F (0)  =  0 

Combine  Equations  (41),  (42),  and  (43)  and  obtain: 
dP 2  /  2  f  2 

-Ji.d  f  - iia dt  +  j-  /  fuwt  +  fe 

«  23  J  U-t|  26S  J  2@ 


(43) 


I 

P  -  /  f(t)dt 

0 _ 

i*  -  i 


(44) 


Term  4>  : 
c 


$  *  sin  a  cos  e 

c  _2- 

p  r 


-.2 


9 

/ 


f(t)dt 


LO 


(45) 


differentiation  gives: 


Si 


2 

=  sin  a  cos  e  f  f  (t)dt 

i  r  J 


(46) 


Now  combine  Equations  HI),  (33),  (37),  (44),  and  (46)  to  obtain: 


♦-(S.^S)  *  62  in  r  +  B2  in -  - - 

X  2/  5(£*  -  e) 


'  I  7  *fe>  dt 
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r 


-  rf*  • 


I 


2|f(S) 


f(t)dt 


2  p  -  /*  f(t)dt 

+  2fh)  _  F7! 


f(5): 


+  sin  a  cos  0  f  f(t)dt 

&  0 


Equate  Equations  (26)  and  (47)  to  font  an  integral  aquation  for 


(i-O 


(r  +  Dm*  (7  +  dm* 

00  00 


2-2 
=  6  in  r  +  5  in 


2/eu*  -  e) 


5  (I*  -2p  f  frrV„  ,  62  f  f(t)  -  f( 

2{f(5)  (5*  -  I)  J  f(t)dt  +  2f(0  J  ,  , 

0  0  1 1  “  t( 


f  WA)  I  f (t)dt 


2  [fj-/ef(t)dt 

+  -£ - 0 

2f(e)  t*  -  t 


On  the  body. 


2  £ 

+  sin  a  cos  0  f  f(t)dt 

*  0 


?  =  R(e) 


Substitute  Equations  (4)  and  (20)  into  Equation  (49)  to  obtain  an 
expression  for  r  on  the  body  in  terms  of  £: 


on  the  body 


9 

|  /  f(t)dt 


Substitute  Equation  (50)  into  Equation  (48),  multiply  through  by  f({), 
and  rearrange  terms  to  obtain: 
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f3(l) 


(l  -  M3)  f(0  +  (7  +  1)M2  -f(|)62^n/|  f  f(t)dt 

V  0 


' £ll)t2  / f(t,dt 


-  ^  <.«-!-/•  £(t)dt 

2  u  -  ti  25  / 


2  /P  -  J5  f(Odt 

6  0 _ 


2  \  (i*  -  i) 


-  2f(|)&  sin  a  cos  6 


Now  factor  a  6  /2  term  from  some  of  the  terms  on  the  right-hand  side  of 
Equation  (51),  combine  the  two  logarithmic  terms,  and  obtain,  finally: 


f3U) 


■  [(■  -  i)  ■ 


25  (7  +  l)Mw  sin  a  cos  e  f(|) 


2  /  6 fi  f(t)dt 

+  ^-(7+  1)M2  -2f  (£)  f  n  (  0  - - 

2  “  \2V  IU*  -  i) 


1(1*  -  t) 


i  i* 

j  f(t)dt  -  f 


5  “  t| 


■\  / £<t)at  ■  ( 


P  -  /  f(t)dt 

0 _ 

(6*  -  I) 


The  boundary  conditions  for  Equation  (52)  are  a3  follows: 

a)  The  cone  apex  (£  =  0)  is  assumed  to  be  a  stagnation  point  for 
zero  and  small  angles  of  attack.  This,  the  total  velocity  (i.e.,  the 


*  *  . 


A. 


sum  o£  the  free-stream  velocity  and  the  perturbation  velocity)  has  to 
be  zero .  Hence , 


U 

00 


+  U  4>-(0)  =  0 

oo  X 


or 


*-(0)  =  -1 


(53) 

(54) 


It  must  be  noted  that  Equation  (54)  Is  only  approximately  true  since  at 
the  stagnation  point  the  small  disturbance  assumption  Is  violated. 

Hence,  the  present  solution  Is  not  correct  at  the  stagnation  point. 
However,  for  a  sharp  slender  cone,  the  stagnation  region  is  very  small 
and  the  error  in  the  vicinity  of  the  stagnation  region  does  not  affect 
the  solution  over  the  rest  of  the  cone.  Substitution  of  Equation  (54) 
into  Equation  (2)  gives  the  boundary  condition  on  f(g)  at  the  stagnation 
point,  namely, 


f(0) 


2-11/2 

+  (7  +  DM.J 


(55) 


b)  The  shoulder  of  the  cone  (g  =  g*)  has  been  assumed  as  the 
location  of  the  sonic  velocity.  It  is  here  that  the  governing  equation 
changes  from  elliptic  to  parabolic  and  hence  it  is  expected  that 

Hi*)  =  0  (56) 

Now  it  Is  observed  that  a  solution  to  the  integral  equation  [Equa¬ 
tion  (52)]  is  dependent  on  a  dimensional  parameter,  g.  Since  very 
little  is  known  about  this  quantity  except  that  it  varies  from  zero  at 
x  =  0  to  g*  at  x  =  1,  it  is  necessary  to  nondimensionalize  Equation  (52) 
with  respect  to  g*  in  order  to  have  an  independent  variable  which  varies 
between  zero  and  one . 

The  following  transformation  nondimensionalizes  Equation  (52)  with 
respect  to  g*  to  make  it  possible  to  solve  for  f(g)  and  at  the  same 
time  eliminates  0  entirely  from  the  equation,  thus  allowing  real 
solutions  for  f(g)  up  to  the  supersonic  shock  attachment  Mach  numbers. 
Let 
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where  £**  is  the  value  of  £/£  at  the  shoulder  of  the  cone  such  that 
0  *  1  £  1.  Thus, 


and  at  the  shoulder 


The  chain  rule  gives: 


t  =  f3£**£ 


§*  =  p§** 


d  _  6_  1  d 

SI  -  a? '  si**  a-e 


d52  J!  dl 


'  ^  IS  («**  as)  ’  /v 


2  -2 
5**  I 


This  transformation  transforms  Equation  (7a)  to: 


(ha«  ffla  <t> 

_li  +  *  +  r  .  ^ee  n 

+  \r  +  “  +  —  =  0  > 

5**  r  r 


and  Equation  (2)  becomes: 

f2(i)  • 


■  l1  -  »!)  -  O’  +  »►£  ♦; 


there fore. 


{l  •  «l)  - 


(7  +  dm; 


Equation  (4)  can  be  rewritten; 


f  (5)  -  i  ~  . 


»  I 


and  integration  of  Equation  (65)  gives: 


3 

x(!)  =|**  J  f(s)ds 


Now  substitute  Equations  (58)  and  (59)  into  Equation  (52)  to  get 
an  integral  equation  which  can  be  solved  for  f(|)  by  numerical  itera¬ 
tion  methods. 


f3(b 


25  h  +  1)M  sin 
00 


a  cos  ej 


„2  , 

+  (7  +  D*r  -2f(l)  ii 


6 J*  f(t)8|**dt 


,2/  0|**5(PS**  -  p|**l)> 


+  f  £(-t)M**  d-t 

P5**t(PS**  -  55**!)  ' 

i  2 

-  f  — dt - —  f  f(t)pg**dt 

q  |p|**5  -  p|**t|  Pft**6  o 


p  -  /5  f(t)p4**  dt 


/P-/ 
\  (st 


(55**  -  P5**l) 


Cancel  5  and  5**  terms  and  the  final  form  of  the  integral  equation  for 
f(5)  becomes: 


f3(!) 


•  [(■  •  o 


25(7  +  1)*^  sin 


a  cos 
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r  /6/"5  f(t)  dt 
+  4  * +  1)M«  [-Mrt)  tnV2°y-s(i~ 


+  Sl-Jh  f  fCod-c-  / 

5(1  -  6)  o  o 


6  -  t 


-  f  f (t)dt  - 
6  n 


g  •  C Hl)di 
(1  -  6) 


Where 


keep  in  mind  that 


at  the  shoulder 


thus, 


■  [(■  •  o 


+  (y  +  Umm 


f(i)  =  o 


3 

x($)  =  £**  j  f(t)dt 


x(0)  =  0 


:(1)  =  | **  f  f(t)dt  =  1 


I**  =  “r - 

J1  f(t)dt 
0 


J. 


A.  ^ _ 


repeating  Equation  (64), 


i-il- 


f2(!) 


(r  +  Dm* 


(75) 


and,  finally,  the  pressure  coefficient  on  the  surface  of  the  cone  is 
given  by  [ 8] : 


C  =  -2*-  -  52  +  (1  -  4  sin2  6)  sin2  a  .  (76) 

P  x 

The  pressure  coefficient  on  the  cone  is  then  determined  as  in  the 
following  paragraphs. 

Equation  (68)  is  solved  for  f(£)  by  numerical  iteration  using  the 
Wu-Aoyama  solution  as  the  first  approximation  The  iteration  is  halted 
when  a  particular  test  condition  is  satisfied. 

The  values  of  f(g)  of  the  last  iteration  are  then  used  to  calculate 

<t>-  according  to  Equation  (75),  and  the  pressure  coefficient,  C  ,  is 
x  p 

calculated  from  Equation  (76). 

Solutions  for  C  can  be  calculated  for  all  Mach  numbers  from  zero 
P 

to  the  shock  attachment  Mach  number. 

A  more  detailed  explanation  of  the  numerical  iteration  technique  is 
presented  in  Appendix  A,  along  with  a  computer  program  listing  in 
Append ix  B . 

The  present  procedure  is  more  consistent  than  that  used  by  Wu  and 
Aoyama  because  both  the  source  strength,  F(|),  and  the  velocity, 
^"(S*^,©),  were  obtained  from  Equation  (11).  Wu  and  Aoyama,  on  the 

other  hand,  used  the  approximate  equation  [Equation  (11)]  to  determine 
the  source  strength,  F(§)>  but  then  used  the  more  exact  equation 
[Equation  (10)]  to  solve  for  the  velocity. 


Section  III.  CALCULATION  OF  FOREBODY  PRESSURE  DRAG 


For  zero  angle  of  attack,  the  pressure,  p,  at  an  arbitrary  section, 

x  (Figure  2),  acts  on  the  projected  area,  dS(x),  to  produce  an  elemental 
force  in  the  axial  direction.  The  free-stream  pressure,  p  ,  acts  at 

the  base  of  the  cone-cylinder  combination  to  produce  a  force  in  the 
opposite  direction.  Hence,  the  total  drag  on  the  forebody  can  be  cal¬ 
culated  by  integration  of  the  net  pressure  (p  -  p^),  which  acts  on  each 

elemental  area,  dS(x). 


,■/ 


(P  -  P  )dS(x) 


x=0 


Now 


dS(x)  =  2*R(x)dR(x)  =  2rtR(x)  — ffl-  dx 

dx 


and 


dfi(x) 

dx 


tan  S  a  6 


for  small  angles;  hence, 


R(x)  =  5x 


(77) 


(78) 


(79) 


(80) 


and 


dS(x)  =  2n&  xdx 

The  net  pressure  force  can  be  expressed  in  terms  of  the  pressure 
coefficient 

p  -  p 

«o 


(81) 


(82) 


f 
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from  which 


P  -  Poo  =  V»  •  <83> 

Therefore,  substitution  of  Equations  (81)  and  (83)  into  Equation  (77) 
gives: 


x=0 


(84) 


Now  the  drag  coefficient  based  on  the  base  area  of  the  cone  is: 


(85) 


where  the  base  area  of  the  cone,  S(i),  can  be  expressed  as 

S(l)  =  JtR2  (1)  =  n&2 


(86) 


Hence,  Equations  (84),  (85),  and  (86)  give: 

1 

Cn  =  2  /  C  xdx 
D  J  p 

0 


(87) 


Liepmann  and  Roshko  [1]  shows  that  for  slender  bodies  of  revolution 
at  small  angles  of  attack,  the  cross  flow  contribution  is  approximately 
2 

equal  to  a  . 

Hence,  the  drag  coefficient  for  small  angle  of  attack  is: 


C 


D 


? 

+  a 


(88) 
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Section  IV.  COMPARISON  OF  PRESENT  METHOD  WITH  EXPERIMENTAL 
RESULTS  AND  OTHER  THEORIES 

Surface  pressure  distributions  calculated  by  the  present  method  for 
a  5-degree  half-angle  cone  are  compared  with  other  theoretical  results 
and  experimental  data  in  Figure  3.  6  is  assumed  equal  to  zero  when  not 

specified.  The  results  labeled  SP1I  were  calculated  by  the  method  due 
to  Smith  and  Pierce[2]  for  a  sharp  cone  with  a  cylindrical  afterbody 
having  the  same  diameter  as  the  base  of  the  cone.  The  results  labeled 
SP1  were  calculated  by  the  same  method  for  a  sharp  cone  with  a  cylindrical 
afterbody  having  a  diameter  approximately  half  as  large  as  the  base  diam¬ 
eter  of  the  cone.  Note  that  the  present  results  are  in  good  agreement 
with  the  experimental  data  for  all  of  the  Mach  numbers,  although  for 
Mach  numbers  near  one  the  calculated  values  are  a  little  lower  than 
the  data  points.  The  previous  Wu-Aoyama  results,  on  the  other  hand, 
are  quite  high  for  the  low  subsonic  Mach  numbers,  but  agree  with  the 
experimental  data  for  the  transonic  Mach  number  range,  except  near  the 
shoulder.  The  error  near  the  shoulder  and  the  high  values  in  the  sub¬ 
sonic  regime  are  probably  due  to  the  fact  that  this  method  neglects  a 
term  in  the  integral  equation  used  to  evaluate  the  pressure  coefficient. 

Consider  now  the  two  subsonic  compressible  calculations.  Note  that 

the  results  labeled  SPI  agree  with  data  for  the  lower  subsonic  Mach 

numbers,  while  the  results  labeled  SPII  agree  better  than  SPI  for  the 

Mach  numbers  above  M  =0.6.  This  is  probably  due  to  the  fact  that  as 

00 

the  Mach  number  Increases,  the  flow  on  the  cone  surface  takes  on  a 
higher  kinetic  energy,  which  helps  the  flow  over  the  shoulder  and 
causes  the  streamlines  to  behave  as  if  the  cone  had  a  full-size  cylin¬ 
drical  afterbody,  regardless  of  the  diameter  of  the  real  afterbody. 

Hence,  one  must  be  careful  when  using  the  subsonic  solutions  because  of 
the  significant  differences  in  pressure  distributions  which  result  for 
different  geometry  inputs.  The  major  disadvantage  of  the  subsonic 
theories  is  that  they  cannot  be  used  for  Mach  numbers  greater  than  or 
equa 1  tc  one . 

Figure  4(a)  shows  the  trend  of  the  SPI  pressure  coefficients  with 
Mach  number,  and  Figure  4(b)  shows  the  same  trend  for  the  SPII  values. 

Note  that  as  the  Mach  number  increases,  the  curves  seem  to  rotate  about  a 
common  axial  location  (x/L  =»  0.72  for  SPI  and  x/L  »  0.9  for  SPII). 

Figure  4(c)  shows  the  behavior  of  the  Wu-Aoyama  solution  with  Mach 
number.  As  Mach  number  increases,  the  curves  simply  translate  to 
higher  pressure  coefficient  values.  Figure  4(d)  shows  the  variation  of 
present  method  with  Mach  number.  As  Mach  number  increases,  the  curves 
seem  to  translate  upwards  as  well  as  rotate  to  a  flatter  slope,  approach¬ 
ing  the  value  calculated  at  the  shock  attachment  Mach  numbers.  From 
Figure  3  the  trend  of  the  experimental  data  is  seen  to  be  similar.  The 
fact  that  the  present  method  exhibits  both  the  rotation  of  the  subsonic 
solutions  and  the  translation  of  the  transonic  solutions  is  very 
encouraging,  and  is  necessary  for  a  unified  theory. 
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Figure  3.  (Continued). 
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Figure  3.  (Continued), 
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Figure  3.  (Continued) . 
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Figure  5  shows  further  comparisons  of  the  present  method  with  data 
and  the  previous  Wu-Aoyama  calculations.  Again,  the  present  solutions 
show  acceptable  agreement  with  data  for  the  subsonic  Mach  numbers; 
however,  the  calculated  values  of  both  the  present  and  the  Wu-Aoyama 
procedures  are  somewhat  lower  than  the  data  points  for  Mach  numbers 
greater  than  one. 

Figure  6  shows  good  agreement  of  the  present  method  with  data  for 
a  15-degree  half-angle  blunted  cone.  Figure  7  demonstrates  Its  applica¬ 
bility  for  the  case  of  a  blunted  cone  for  a  6-degree  angle  of  attack. 
Figure  8  Is  further  proof  that  the  method  gives  good  results  at  moderate 
Mach  numbers  for  small  angles  of  attack. 

It  Is  Interesting  to  note  the  variation  of  the  pressure  coefficient 
with  Mach  number  at  certain  fixed  axial  locations  on  the  cone  surface. 
Figure  9  points  out  the  advantages  the  present  method  has  over  the 
other  calculated  solutions.  A  comparison  of  the  solutions  for  three 
different  positions  is  shown  in  Figure  10. 

Also  of  interest  is  a  comparison  of  the  present  solution  with  other 
transonic  solutions,  as  shown  in  Figure  11.  It  must  be  noted  that  only 
the  present  and  Wu-Aoyama  solutions  were  specifically  intended  for 
Mach  numbers  other  than  one.  The  present  method  exhibits  the  best 
shape,  although  it  seems  to  fall  below  the  data  points. 

A  drag  coefficient  versus  Mach  number  curve  is  shown  in  Figure  12 
for  two  different  cone  angles  . 
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Figure  5.  Comparison  of  present  method 
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Figure  5.  (Continued). 
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Figure  8.  Surface  pressure  distribution  for  s  sharp  cone  at 
a  4-degree  angle  of  attack. 
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Variation  of  C  with  free-stream  Mach  number  for  a  sharp  cone 
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Figure  12.  Comparison  of  pressure  drag  with  Mach  number  (present  method). 


Section  V.  CONCLUSIONS  AND  RECOMMENDATIONS 


It  can  be  concluded  from  this  study  that  the  present  method  provides 
the  first  real  unified  solution  to  the  flow  on  the  surface  of  a  cone. 

It  is  applicable  for  the  full  range  of  Mach  numbers  right  up  to  shock 
attachment.  However,  because  the  shape  of  the  sonic  line  is  not  known 
a  priori  as  a  function  of  Mach  number  and  cone  angle,  the  total  flow 
field  solution  cannot  be  computed.  Although  the  solutions  exhibit  the 
correct  tendencies  with  Mach  number,  the  values  calculated  above  M  = 

0.95  are  somewhat  lower  than  the  experimental  results.  This  is  due  to 
the  fact  that  the  solution  is  only  a  first  approximation.  One  recom¬ 
mendation  would  be  to  apply  a  successive  approximation  technique  to 
solve  for  the  complete_velocity  potential.  This  deviation  may  also  be 
a  result  of  forcing  f(|)  positive  when  solving  the  integral  equation 
[Equation  (68)].  Another  recommendation  would  be  to  allow  f(£)  to 
assume  negative  values,  iterate  until  the  solution  converged,  even 
though  the  maximum  value  of  x(§)  would  be  greater  than  one,  and  then 
compress  the  x(§)  scale  about  the  center  x(f)  such  that  0  £  5t(g)  £  1. 
Preliminary  calculations  indicate  this  would  probably  raise  the  pressure 
coefficient  values  all  along  the  cone  surface. 

Different  rearrangements  of  the  terms  of  Equation  (68),  as  dis¬ 
cussed  in  Appendix  A,  might  also  be  explored. 

Consider  also  the  possibility  that  viscous  effects  become  important 
for  these  Mach  numbers.  Another  recommendation  would  be  to  correct  the 
present  inviscid  solutions  for  viscous  effects  in  an  effort  to  find 
better  agreement  with  experimental  data.  Finally,  it  is  recommended 
that  the  shape  of  the  sonic  line  be  studied  in  an  attempt  to  obtain 
solutions  for  the  whole  flow  field. 

Although  all  of  these  recommendations  relate  to  adjustments  of 
the  present  method,  the  possibility  that  the  experimental  data  are 
inaccurate  for  those  upper  transonic  Mach  numbers  should  not  he  com¬ 
pletely  ruled  out  either.  It  is  a  well-known  fact  that  wind-tunnel 
wall  interference  effects,  as  well  as  compressibility  effects,  become 
very  significant  in  this  Mach  number  range. 
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Appendix  A.  DETAILS  OF  THE  NUMERICAL 
ITERATION  PROCEDURE 


Details  of  the  numerical  procedures  used  in  the  present  analysis 
are  described  below. 

Rewrite  Equation  (68)  as: 

2 

f3(|)  =  A  +  §r-  (7  +  1)M2  (-B  +  C  -  D  -  E  -  F]  (A-l) 

L  oo 


5  .  . 

and  remember  x(g)  =  g**  J  f(t)dt. 

0 

The  values  of  f(l),  g**,  and  x(g)  calculated  by  the  original  Wu- 
Aoyama  method  are  used  as  the  first  approximation  values  in  the  present 
iteration  scheme,  and  are  substituted  into  Equation  (A-l)  to  solve  for 
•• 

new  values  of  f(g).  A  new  value  of  g**  is  then  calculated  according 


to  4**  =  1//*  f(l)d|  using  the  new  values  of  f(4),  and,  finally,  new 
0 

values  of  x(|)  are  calculated  according  to  x(g)  =  £**  f(|)d|. 

0 

Several  tests  are  then  made  to  see  if  the  new  values  of  f(g)  are 
acceptable  and  satisfy  convergence.  If  they  do  not,  then  these  new 
values  of  §**,  f(f),  and  x(|)  are  substituted  into  Equation  (A-l)  to 
begin  a  second  iteration.  Successive  iterations  are  performed  >mtil 
the  solution  converges. 

A  plot  showing  the  relative  magnitudes  of  the  individual  terms  of 
Equation  (A-l)  is  presented  in  Figure  A-l.  The  terms  were  calculated 
during  the  first  iteration  using  the  Wu-Aoyama  values  of  x(|),g**,  and 
f (4) .  It  is  easy  to  see  that  term  B  is  dominant  except  near  the  shoulder 
(4  =  1)  where  term  C  becomes  dominant.  Term  A  is  not  shown  because  it 
is  small  for  the  particular  case  being  discussed  and  varies  from  0.0175 
at  4  *  0 .01  to  0 .0113  at  4  =  0.99.  At  lower  Mach  numbers  the  value  of 
term  A  would  increase  somewhat. 

It  should  also  be  noted  that  the  factor 


c  2  2 

\  (y  +  DM, 


is  very  dependent  on  the  cone  angle  and  the  Mach  number.  Even  sj, 
very  close  to  the  shoulder,  the  right-hand  side  of  Equation  (A-l) 
becomes  negative.  This  is  caused  by  term  C  which  has  a  tendency  to 
blow  up  at  4  =  1.  This  singularity  is  artificially  eliminated  in  the 
computer  program,  as  follows. 

As  long  as  the  right-hand  side  of  Equation  (A-l)  is  positive,  the 
value  of  f(^)  is  calculated  according  to  Equation  (A-l).  However,  when 
the  right-hand  side  becomes  negative,  the  value  of  each  successive  f(g) 
is  forced  to  equal  one-half  the  value  of  f(4)  of  the  previous  iteration. 
For  example,  if  f(0.89)  ■  0.1  and  f(0.90)  »  -0.2,  as  calculated  by 
Equation  (A-l),  then  f(0.90)  would  be  set  equal  to 


0.5  f (0.90) 


previous  iteration 


f (0.91)  =  0.5  f (0.91) 


previous  iteration 


etc.,  up  to  the  shoulder.  In  effect,  this  procedure  is  simply  limiting 
the  effect  that  term  C  can  have  on  the  solution. 


-  A 


Figure  A-l.  Comparison  of  terras  of  integral  Equation  (A-l) 


/} 


Now  there  are  two  tests  for  convergence: 

1)  Test  No.  1:  Whenever  the  maximum  absolute  difference  between 
the  new  value  of  f(g)  and  the  previous  value  of  f(g)  is  less  than  or 
equal  to  0.005,  the  iteration  is  halted.  This  works  for  the  low  Mach 
number  cases.  However,  for  Mach  numbers  in  the  transonic  regime,  it 
was  found  that  the  Iterative  solution  developed  an  instability  near  the 
shoulder  if  the  number  of  iterations  became  too  large. 

2)  Test  No.  2:  Hence,  the  iteration  was  stopped  as  soon  as  the 
pressure  coefficient  values  achieved  a  monotonically  decreasing  trend 
up_to  a  value  of  x/L  =  0.99.  Figures  A-2  and  A-3  show  typical  plots  of 
f(g)  as  calculated  by  the  present  method.  A  minimum  of  three  iterations 
was  always  calculated  after  the  Wu-Aoyama  solution  was  obtained. 

A  solution  was  said  to  be  nonconvergent  if,  after  10  iterations  of 
the  Wu-Aoyama  preliminary  calculations,  test  No.  1  had  not  been  satisfied, 
or  if  after  25  iterations  after  the  Wu-Aoyama  solution  was  obtained,  both 
tests  No.  1  and  No.  2  were  still  not  satisfied. 

Convergence  was  found  to  be  very  dependent  on  Mach  number. 

Figure  A-4  shows  the  maximum  Mach  numbers  for  which  solutions  converged. 


Now  the  velocity,  $-,  is  calculated  according  to  Equation  (75)  and, 

hence,  values  for  the  pressure  coefficient,  C^,  are  obtained  using 

Equation  (76).  Then  the  local  Mach  numbers  at  each  location  on  the 
surface  of  the  cone  are  calculated  according  to: 

r  ii/2 


(A-2) 


which  was  derived  using  the  isentropic  flow  conditions  where  the  total 
to  static  pressure  ratio  has  been  expressed  in  terms  of  the  pressure 
coefficient  as  follows: 


and 
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CONE  SEMI-APEX  ANGLE,  6  |d»«) 

Figure  A-4.  Comparison  of  shock  attachment  Mach  numbers  of 
present  method  with  supersonic  flow  theory. 


Also,  the  free-stream  to  static  pressure  ratio  is  calculated  for 
all  points  on  the  cone  surface  according  to: 


r00 

p 


C  r  M  +  1 

p  2  « 


Finally,  the  pressure  drag  coefficient  is  calculated  for  the  zero 
angle  of  attack  cases  according  to: 


J 


C 

P 


set)  ^ 

di 


where 


— =  i**  f(i) 

ai 


It  should  be  noted  that  several  different  techniques  were  tried  in 
order  to  cause  the  f(i)  values  to  drop  towards  zero  closer  to  the 
shoulder,  as  shown  by  the  dotted  line  in  Figure  A-5.  However,  since 
x(i) ,  £**,  and  f(i)  are  all  dependent  on  one  another,  the  only  effect 
was  that  £**  assumed  a  different  value  for  the  new  f(i)  values  and 
left  the  surface  pressure  distribution  relatively  unchanged. 

Another  trial  included  letting  f(|)  assume  negative  values.  This 
resulted  in  values  of  x  greater  than  one.  Although  this  technique  was 
not  pursued  in  this  study,  it  was  noted  that  a  more  thorough  investiga¬ 
tion  may  prove  worthwhile  . 

The  x/L  values  can  be  collapsed  onto  a  compressed  scale  such  that 
0  ^  x/L  ^  1,  which  might  adjust  the  surface  pressure  coefficient  values 
closer  to  those  obtained  by  experiment. 

Another  technique  that  might  be  used  is  the  rearrangement  of  the 
terms  on  the  right-hand  side  of  Equation  (68)  to  see  if  the  effect  of 
the  singularity  at  the  shoulder  can  be  minimized. 
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%,  v,  0 


Figure  A-5.  Typical  plot  of  f(g)  near  the  shock  attachment 
Mach  number. 
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Appendix  B.  COMPUTER  PROGRAM  AND  SAMPLE  OUTPUT 


This  appendix  presents  a  listing  of  the  computer  program  which 
calculates  the  surface  pressure  distribution  on  a  cone  according  to  the 
present  method,  followed  by  a  sample  output.  The  program  was  written 
for  the  IBM  370/155  computer  at  Arnold  Engineering  Development  Center, 
Arnold  Air  Force  Station,  Tennessee.  The  drag  coefficient  is  only 
printed  for  the  case  of  zero  angle  of  attack.  To  calculate  the  drag 

2 

coefficient  for  the  case  of  small  angles  of  attack,  simply  add  a  to 
the  zero  angle  of  attack  value  [see  Equation  (88)]. 

A  list  of  symbols  used  in  the  output  along  with  their  meanings  is 
presented  in  Table  B-l. 
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TABLE  B- 1 .  COMPUTER  PROGRAM  SYMBOLS 


Notation  of 
Computer  Output 

Notation  Used  in  Report 

MACH  INF 

M 

00 

CONE  SEMIVERTEX  ANGLE 

6 

GAMMA 

y 

ANGLE  OF  ATTACK 

a 

AZIMUTH  ANGLE 

9 

ITERATION 

Total  number  of  iterations 

performed  (including  those  of 
the  Wu-Aoyama  calculations) 

KSAI  (SONIC) 

5** 

DRAG  COEFFICIENT 

cD 

D1 

KSAI 

1 

X/L 

x(!) 

F 

f(b 

DF 

f^^last  ”  previous 

•iteration  1  iteration 

U 

X 

CP 

C 

P 

M 

M 

X/D 

Ratio  of  axial  position  on  cone 
to  base  diameter  of  cone 

|  PI/P 

P« 

P 
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CO*MnuC 
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I 

i 
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N 
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! 

•- 

\ 

i 

! 
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fcl 

t 

X 

2 

■* 


w 

I 


I 

I 
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t 
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SYMBOLS 


A,B,C,D,E,F 


Terms  of  Equation  (A-l)  as  described  in  Appendix  A 
Pressure  drag  coefficient 

Pressure  coefficient 

Pressure  drag 
Source  strength 

Function  defined  by  Equation  (43) 
transformation  function  defined  by  Equation  (2) 
Function  defined  by  Equation  (12) 

Function  defined  by  Equation  (40) 

Function  defined  by  Equation  (5b) 

Length  of  cone 
Mach  number 

Static  pressure  on  cone  surface 
Total  pressure 


Dynamic  pressure 
Body  radius 
Radial  coordinate 
Projected  area 
Dummy  variable 
Dummy  variable 
Free-stream  velocity 


Perturbation  velocity  in  axial  direction,  <t>-  U 

X  co 

Perturbation  velocity  in  radial  direction,  $-  U 

r  oo 

Perturbation  velocity , (®  /r)U 

Z7  00 

Axial  coordinate 
Angle  of  attack 
/ - Z 

y  1  -  M  =  compressibility  factor 
Ratio  of  specific  heats 
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€ 

t 

9 


Subscripts 

a 

c 

1 

r 

t 

x 

0,1, 2, 3 

5 

e 

oo 


Superscripts 

(') 

* 

** 

(  )' 


Half  angle  of  cone 

Small  number  defined  in  Equation  (5b) 
Independent  variable  defined  by  Equation  (4) 
Roll  angle 

Linearized  perturbation  velocity  potential 


Axial  flow 
Cross  flow 

Axial  flow  (drag  calculations) 
Indicates  differentiation  with  r 
Total  condition 

Indicates  differentiation  with  x 
Identification  numbers 
Indicates  differentiation  with  § 
Indicates  differentiation  with  3 
Free-stream  condition 


Indicates  nondimens ionali zed  quantity 
Conditions  at  the  sonic  point  on  the  cone  surface 
Conditions  at  the  sonic  point  on  the  cone  surface 
Differentia  tion 
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